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Abstract 

In this study, the particles of the quantum gases, namely bosons and 
fermions are regarded as g-ons by the parameter of the fractional exclu- 
sion statistics g, where < g < 1. With this point of departure, the 
distribution function of the g-on gas is obtained by the variational, steep- 
est descent and statistical methods. The distribution functions which are 
found by means of these three methods are compared. It is seen that the 
thermostatistical formulations of quantum gases can be unified. By suit- 
able choices of g standard relations of statistical mechanics of the Bose 
and Fermi systems are recovered. 



Introduction 

Particles of quantum systems, including well known fermions and bosons 
can be called as g-ons which obey fractional exclusion statistics. The aim of 
this study is to devolop the distribution function of a g-on gas by means of 
variational, steepest descent and statistical methods and discuss them. Futher- 
more, another purpose is to unify the different thermostatistical formulations 
of quantum gases where the concomitant distribution function of g-on gas is 
used. 

The outline of this paper is as follows: In section 1.1, Haldane statistics 
is obtained by the variational method where the entropy of the g-on gas is 
employed. By means of an approximate solution to Haldane statistics, an elegant 
formula for the distribution function of the g-on gas is found. In section 1.2, 
where the steepest descent method is applied to the partition function, an exact 
solution of distribution function is obtained. In section 1.3, the probability of 
occupation numbers and average distribution functions are investigated by the 
statistical method. In section 2 and 3, fluctuations in the occupation number 
and thermodynamics of the g-on gas are respectively discussed. 
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1 Derivation of the distribution function of the 
g-on gas system 



1.1 Variational method 

Let us consider an ideal g-on gas system which consists of N particles. Let 
the system be in heat and particle baths that evolves towards equilibrium, that 
is, the total energy E and the total number of g-ons N fluctuate, but, these 
quantities are conserved on the average. The entropy of the ideal g-on system 
which is not in equilibrium has been obtained from the statistical weight and is 
given by [1], 



K 

S g = 9k [1 + n k (1 - g)\ ln[l + n k (1 - g) - (1 - gn k ) In [1 - gn k ] - n k \nn k 
fe=i 

(1) 

where the Boltzmann constant is set to unity, g k is the degeneracy, n k is the 
occupation number, K is the total number of states; g is the parameter that 
characterizes the complete or partial action of the exclusion principle i.e. makes 
interpolation between quantum gases [2,3]. Following the variational method 
outlined in the ref.[l], the following equation is obtained 

[In 1 + n k (1 - g)f- 9) + (1 - g) + In (1 - gn k ) 9 + g - lnn fe - 1 - a - f3e k = 

(2) 

where a = — £ , (3 — ^ , e k is is the energy of the k. state and pi is the chemical 

potential. After performing some algebra Eq.(2) becomes; 

[l + nai-g)] 1 ' 9 {l-gn k f = &Xk 
n k 

where 



x k = (3 (e k - n) . (4) 

Eq.(3) is known as Haldane statistics [4]. 

Let us write Eq.(3) in a more compact form 

- (1 - V) 9 = e x (5) 

y 

where y is defined as 

y = ^~i — • (6) 
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In order to get a real x on the right hand side of Eq.(5), y < 1. This inequality 
imposes a condition 

gn k < 1. (7) 

Note that for g=l inequality(7) gives rife < 1 , for g=0; can take any val- 
ues and for the intermediate values rife is restricted asm < K Therefore, in- 
equality^) could be interpreted as the generalization of the Pauli Exclusion 
principle. 

Since it is difficult to solve Eq.(5), an approximation is needed for y . Ex- 
panding the factor (1 — y) 9 in a Taylor series about y =0: 

(1 - y) 9 nl-gy + O (y 2 ) . (8) 

Thus, Eq.(5) becomes 

nk = e /3( £fc -M)+2.g-r (9) 

The last approximation is a linear approximation to Eq.(3) and therefore does 
not satisfy inequality (7), if g is not close to or 1. Interesting enough, Eq.(9) 
turns out to be the Maxwell distribution when g= \ [1] ■ 



1.2 Steepest Descent Method 

In order to get another solution of distribution function nk we are willing to 
calculate it using the steepest descent optimization method [5]. Two of the 
authors of the present study FB and DD recently introduced a partition function 
Z g for a grandcanonical ensemble of a g-on gas [1]. Let us write it for the 
extensive systems as 

K ( 1 oo ^ 

2 9 =n E^+i 1 -^^ 1 ( io ) 

fc=l tfe=0 fe=2 J 

where Xk is given by Eq. (4) . The given partition function leads to the following 
partition function for the canonical ensemble 



K 

z 9 = II (! + + (! - a) {vt + vl + ••■) (ii) 

k=l 

where yt = e^ l3ek and k denoting the quantum state £fe. Here g is a parameter 
that unifies the partition functions of boson and fermion systems under g-on 
gas system. For g=0 the partition function in Eq.(ll) becomes the partition 
function of a boson system, for g=l it becomes that of a fermion system. Let 
us find the distribution function of a g-on system using this partition function 
and the steepest descent method. 
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In order to use the steepest descent method one has to express the partition 
function as a complex integral around a closed contour. For this purpose let us 
define the function 

h (*) = ft {(1 + *Vk) + (1 - .9) (zVk + zhl + -..) } • (12) 

k=i 

Using this function, the partition function given in Eq.(ll) can be written 

as 

c 

or 

*. / n + — zVl + ~ yl + - ) )] »■ <«) 

c 

This integral (or partition function) can be approximated by using the method 
of steepest descent. In order to see it let us write 



e 



h(z) 



(15) 



where f g (z) is given by Eq.(12). Using this definition one can rewrite the h(z) 
as 



K 



h(z) = lnll 



k=l 



1 + zx k (1 - g) z 2 y 2 k (1 + + z 2 y 2 k 



yN+l 



(16) 



or as 



h(z) = {In [(1 + zy k ) + (l-g) z 2 y 2 (l + zy k + z 2 y 2 + ...)] } - lns"^ 1 . 

(17) 



fc=i 



The function h(z) has a minimum at z and this minimum is very steep 
because the function f g (z) consists of infinite number of multiplication and 
N + 1 is very large. Thus by using the approximation formula of the steepest 
descent method, one gets an expression for the partition function of a g-on gas 
as [6] 
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Z g =[zf +1 f g (z Q )2irf"(z )] 1/2 . (18) 

In order to find the distribution function one applies the steepest descent con- 
dition to h(z), that is 

h, (Zq) = h > j- r — oo/t 9~^5 7 = U - 

z ° i~o 1 + Z ° Vk + '4v k (! + z oVk + z&vt + •••) 

(19) 

Since N » 1, one can write N+l « N. Using this approximation and arranging 
Eq.(18) one obtains 

N „ y> zoy fc + (l-. 9 ) (2z 2 ^ + 3z 3 ^ + ...) 

~ fro 1 + Zoyfe + ~ ^ ^ + z ° Vk + z o2/l + -) ' 

The factor can be interpreted as fugacity e"^ [6]. For the sake of simplicity 
let us define a new variable as 

a k = e-^-^ = e- x * = z y k . (21) 

Using the series 

l + a k +a 2 k + al + ... = (l-a k y 1 (22) 

l + 2a k + 3al + ... = {l-a k )- 2 (23) 

where a k < 1. By means of these series and the definition (21) one can write 
Eq.(20) as 

N = a k = ^ a fc (1 - ggofc + gq|) 



or 



(1 - e~ x f) (1 -ge~ 2x ") 
Since the total number of particles of a system is given as 



N = J2nk-, (26) 



fc=i 



the terms in Eq.(25) gives the number of particles which occupy the k th state. 
Hence for the distribution of g-on particles it is possible to write 

_ e~ Xk (l - 2ge- Xk +ge- 2xk ) 

nk = (l-e- Xk )(l-ge- 2xk ) ' (2?) 
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Eq.(27) which unifies the distribution of the quantum gases recovers the well 
known Bose and Fermi distributions for g=0 and g=l respectively. 

In order to find a more elegant result which could be used moreeasily in the 
calculations the term in Eq.(24) can be written as 

g 4-2a, + l 
1 - 2ga k + ga% 

By taking into account Eq.(28), Eq.(27) reads, 

^ = l + a t °(2 g -l) (M) 
which is valid about g=l and g=0. Let us write Eq.(29) explicitly; 

~ k= e^ k -^ +2g -l (30) 

It is observed that Eq.(30) is same as Eq.(9) which is obtained through an 
approximation to Haldane formula. 

1.3 Statistical Method 

The partition function of g-ons which is given by Eq.(lO) can be used to find 
the probability of the k th state being occupied by n k g-ons. If the state k with 
an energy is occupied by n k g-ons then the other states are occupied by 
N-nfe g-ons, where N is the total particle number. Therefore the ratio of the 
partition function for the states except k th state times e~ nkXk to the partition 
function gives us the probability of finding n k particles in the k. state. In the 
thermodynamic limit i.e. as N — > oo this probability can be written as 



f 1 oo 

e -n kXk Yl J J2 e~ niXi + (1 - g) e- 

p " K) = fi = — ■ ( 31 ) 

H \ E e ~ ntXt + (! - S) E e-"'*' 

i=l l«t=0 n k =2 



or writing in a simplified form: 

p -n k x k 

P 9 k (»*) - = • (32) 

e -n k Xi _|_ [I _ ^2 e~ nkXi 
n k =0 n k =2 

After performing the summation Eq.(32) turns out to be 



= ,_„-L k ■ ( 33 ) 
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It could be mentioned that Eq.(33) is a unified form of the probabilities of 
occupation number of the quantum gases simply bosons and fermions. One 
may write the distribution functions by means of Eq.(33) for g-ons as 



n^= J2 nkP a + t 1 _ -9) E n k P g ( n *) ■ 

n k =0 n k =2 

After performing the summations it reads 



(1 



1 — ge 2xk 
or more explicitly; 

nk 



9 + 



1-9 



(1-e- 



e - x "(ge- 2xk - 2ge~ Xk + 1) 
(1 - ge~ 2xk ) (1 - e- x ") 



e -P(e k -n) _ 2ge~ l3 ( £k ~^ + ge~ 2/3( - Ek ~^) 



(l - e-0( e >>-ri) (l - ge- 2 l 3 ( ek -^) 
which is the same as Eq.(27) as expected. 



(34) 



(35) 



(36) 



2 Fluctuations in the occupation number 

In the preceding section we have calculated the probability Pg (rife) of occupying 
a g-on of the k th state in the Fock space and the distribution function nk . 

The deviations of mean values in the distribution function are given by the 
standard deviations of distributions, 

(37) 



9 o 9 

cr = ni - rife . 



Therefore at first n\ has to be calculated. For this purpose let us write the 
following equation 



1 oo 

A= E n l p " ("*) + a - .9) E n l p " (»*) • 

nfc=0 nfc=2 



(38) 



Substituting the expression for P k (n k ) into Eq.(38) and performing the sum- 
mations one finds 



— e~ Xk (1 - e~ x >>) 



1 - ge 



-2x k 



9+ (I- 9) 



1 + e" 



(1-e- 



(39) 



By means of Eqs.(35),(37) and (38) and one can calculates the variance as 



%k ( „2„-ix k 



Age 



-3x k 



+ 6ge~ 2xk -Age~ Xk + 1) 



(1 - e~ x ") (l-ge- 2x "Y 



(40) 



and the ratio =^ as 



a 2 ^ gXfc^g-^fc _4 ge -3xk + Q ge -2x k -4ge~ Xk + 1) 
nfe 2 (c/e" 2 ^ - 2#e" x * + l) 2 

which reduces to the standard results. 
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3 Thermodynamics of the g-on gas 

In this section, as a concrete application of g-on statistics we would like to 
investigate the thermodynamics of an indistinguishable g-on gas. The ideal g- 
on gas is a model system, where the influence of quantum statistical effects can 
be very well studied. Our aim is to calculate the thermodynamical potential of 
the g-on gas $ S (T, V, jj). 



$ 9 (T,V,fj,) = -k B TlnZ g (T,V,z) 



(42) 



where summation approximately is over all of the quantum states. By taking 
into account equation (10) one derives approximately 



K 



$ 9 (T,V,fi) = -k B Tj2^ 



k=\ 



1 



1 



where 



+ .9-1 



(43) 



(44) 



Since it is difficult to handle with exact formula for n k its approximation Eq.(10) 
or Eq.(30) is appropriate 



1 



n k = 

u k + g 

Given the generalized thermodynamic potential 

$ 9 (T,V,») = -p g V 
the pressure of the g-on gas is obtained as 

K 

P g V = k B T In 



k=\ 



1 



while the total number of particles is found as 



E— ■ 

u Uk+9 



(45) 



(46) 



(47) 



(48) 



The state of the gas is determined by Eqs. (47) and (48). But one can also 
calculate the energy of the gas which is given by 



E„ 



u k + g 



The logarithm of the grandpartition function Q g (T,V,z) is given as 

1 

n I + 

fe=i 



K 



Q g (T,V,z) = lnZ g (T,V,z) = J2^ 



(49) 



(50) 
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Let one of the g-on energies Sk are those of the energies of a free quantum 
mechanical particle in a box of volume V. The chemical potential or fugacity is 
not fixed, but the total number of g-ons are conserved on the average and z has 
to be determined from Eq.(48). For a large volume the sum over all single g-on 
states can be given as the single g-on density states, where spin is not taken 
into account. Hence, 

D (e) = Ae 1 / 2 (51) 

can be written where 

A =™{2 m hf 2 . (52) 
Thus the summation in Eq.(50) becomes 

1 



Q g (T,V,z) = J In 




1 + 



z 1/2 de (53) 



where Eq.(51) is taken into account. When the integration is performed by 
parts, Q(T,V,z) is reduced to: 

00 

Q g (T, V, z) = ^A(3 J e^ 2 nde. (54) 


After substituting n from Eq.(45) and then rearranging 



2 f s 3/2 
3 Af3Z J e^-(\-2g)z 



Q g (T, V, z) = -Apz I — p - — ds (55) 







is obtained, where z is the usual fugacity i.e. z — e^ M . The integral is calculated 
using the related formula pp. 326 of reference[7]. Thus one concludes that for 
a g-on system 



Q g (T,V,z) = ^Af3zF 5/2 (g,z) (56) 



where the definition 



Fiqz)= m Y [z{1 - 29)]k (57) 

is adopted in which T (n) is the usual gamma function. Then the state of the 
g-on system is determined by the following equation: 

Pg V = k B TQ g (T, V, z) . (58) 



9 



Thus, the equation of state of fermions and bosons is unified 

In thermodynamics limit, that is in the limit of infinite volume with particle 
density held fixed, particles in the ground state, having no contribution of the 
kinetic energy to the pressure, then, in terms of the thermal wavelength, the 
pressure of the Bose system becomes, as it is expected, 

kT ~ z k 

?o = ^E^ (^) 
fe=i 

where T (|) = 3-^ is taken and Ais the thermal wavelength [8]. 

In a similar manner,the thermodynamical properties of a Fermi system fol- 
lows immediately from the logarithm of the grandpartition function of a g-on 
system, which is given by Eq. (43). In terms of the thermal wavelength, the 
pressure of the Fermi system is 



pi = a. 



k B T ~ (-l) k - 1 z k , s 

T^HvT- (60) 



fe=i 



where g s is a factor due to spin [8]. 

For a large volume, in accordance with equation (49), the total number of 
g-ons could be written as 

°° l/ 2 

N a {T,V,z) = Azf efle _ e * (61) 

o 

assuming that single number of g-on states can be calculated in terms of in- 
tegrals. Calculation of the integral in Eq.(61) using the related integrals in 
reference[7] pp326 leads to 

N g (T,V,z) = AzF 3/2 (g,z) (62) 

where F 3 / 2 (g, z) is defined by Eq.(58) [8]. If g=0 and g=l arc substituted in 
Eq.(62) the formula for the total number of Bose system No (T, V, z) and Fermi 
system N\ (T, V, z) are respectively recovered . 

For a Bose system the total number of particles iV (T, V, z) including e = 
energy states, is found to be: 

N (T,V,z) = AzF 3/2 (z) + (63) 

where the last term represents the contribution of the energy level e — to the 
mean particle number. 

In a similar manner, for the total number of fermions one gets 

-IT OO / .I N k— 1 k 

Nl (T,V,z) = 9s ^J2 fci/2 " ( 64 ) 
fe=i 
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where V (3/2) = \/n/2 is taken. This is an expected standard result [8] 

Conclusions 

This study has been initiated from the point of view bosons and fermions, 
which obey quantum statistics, could be regarded as g-ons which obey fractional 
statistics. In this paper, Haldane statistics is rederived by variational method 
where entropy of the g-on gas is involved. Due to the importance of the dis- 
tribution function, by means of steepest descent and statistical methods, mean 
occupation number is derived where partition function of g-on gas is used. It 
is concluded that mean occupation number which is found by these methods 
lead to same expressions. Within this context the probabilities of occupation 
number of g-on gas is also calculated which unifies the probabilities of quantum 
gases. 

The approximate solutions of the distribution function of a g-on gas system, 
which is found from the variational, steepest descent and statistical methods 
are developed to an elegant expression that recovers FD distribution for g=l, 
and BE distribution for g— 0. It should be remarked that the partition function 
which is used in this study is appropriate for the g-on gas since it leads to the 
same approximate expression for the mean occupation number with Haldane 
statistics. 

As an application of the mean occupation number, the thermodynamical 
formulation of the g-on gas is derived which gives the opportunity to unify the 
thermostatistical formulation of quantum gases, simply fermions and bosons. 

Authors would like to thank Ege University Research Fund for their partial 
support under the Project Number 2000 Fen 19. 
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